A posteriori error analysis of a quadratic finite volume method for nonlinear elliptic problems

In this article, we construct and analyze a residual-based a posteriori error estimator for a quadratic finite volume method (FVM) for solving nonlinear elliptic partial differential equations with homogeneous Dirichlet boundary conditions. We shall prove that the a posteriori error estimator yields the global upper and local lower bounds for the norm error of the FVM. So that the a posteriori error estimator is equivalent to the true error in a certain sense. Numerical experiments are performed to illustrate the theoretical results.     

Second moment of error term of mean value for binary Egyptian fractions

We study the mean square of the error term of the mean value for binary Egyptian fractions.We get an asymptotic formula under the Riemann Hypothesis.     

More efficient time integration for Fourier pseudospectral DNS of incompressible turbulence

Time integration of Fourier pseudospectral DNS is usually performed using the classical fourth-order accurate Runge-Kutta method or other second- or third-order methods, with a fixed step size. We investigate the use of higher-order Runge-Kutta pairs and automatic step size control based on local error estimation. We find that the fifth-order accurate Runge-Kutta pair of Bogacki and Shampine gives much greater accuracy at a significantly reduced computational cost. Specifically, we demonstrate speedups of 2× to 10× for the same accuracy. Numerical tests (including the Taylor-Green vortex, Rayleigh-Taylor instability, and homogeneous isotropic turbulence) confirm the reliability and efficiency of the method. We also show that adaptive time stepping provides a significant computational advantage for some problems (like the development of a Rayleigh-Taylor instability) without compromising accuracy.     

一种非正交对称条带束流位置监测器设计

介绍了一种条带束流位置监测器(BPM)的设计与仿真方法。在国家同步辐射实验室"太赫兹近场高通量材料物性测试系统"工程项目中,针对波荡器出口处真空室非正交对称性的问题,设计了矩形真空室和跑道形真空室下的两种非正交对称性条带BPM,并与传统的圆形真空室下条带BPM进行对比。基于边界元法,利用MATLAB软件分别对三种真空室下的条带BPM进行建模和仿真。仿真结果表明:相对于传统的圆形真空室下条带BPM,矩形和跑道形真空室下条带BPM灵敏度提高了30%,阻抗匹配误差相对降低了20%,束流位置拟合误差降低了80%。考虑加工精度,矩形真空室下的条带BPM更适用于该工程。     

Response surface methodology toward the optimization of high-energy carotenoid extraction from Aristeus antennatus shrimp

High-energy assisted extraction techniques, like ultrasound assisted extraction (UAE) and microwave assisted extraction (MAE), are widely applied over the last years for the recovery of bioactive compounds such as carotenoids, antioxidants and phenols from foods, animals and herbal natural sources. Especially for the case of xanthophylls, the main carotenoid group of crustaceans, they can be extracted in a rapid and quantitative way with the use of UAE and MAE.     

A study of properties of adaptive quadratic grids for three-dimensional near-surface field computations

Curved geometries and the corresponding near-surface fields typically require a large number of linear computational elements. High-order numerical solvers have been primarily used with low-order meshes. There is a need for curved, high-order computational elements. Typical near-surface meshes consist of hexahedral and/or prismatic elements. The present work studies the employment of quadratic meshes that are relatively coarse for field simulations. Directionally quadratic high-order elements are proposed for the near-surface field regions. The quadratic meshes are compared with the conventional low-order ones in terms of accuracy and efficiency. The cases considered include closed surface volume calculations, as well as computation of gradients of several analytic fields. A special method of adaptive local quadratic meshes is proposed and evaluated. Truncation error analysis for quadratic grids yields comparison with the conventional linear hexahedral/prismatic meshes, which are subject to typical distortions such as stretching, skewness, and torsion.     

Terpolymerization of Triisopropylsilyl Acrylate,Methyl Methacrylate,and Butyl Acrylate: Reactivity Ratio Estimation

Ternary monomer reactivity ratios of triisopropylsilyl acrylate (SiA), methyl methacrylate (MMA), and n‐butyl acrylate (BA), as common monomers in self‐polishing coatings (SPCs) binders are obtained using experimental data collected from free radical bulk polymerization at 70 °C. Different terpolymerizations at low and medium‐high conversions are performed at optimized feed compositions. Estimations are made using the error‐in‐variables model (EVM) framework, applying the recast form of the Alfrey–Goldfinger (AG) model and a direct numerical integration (DNI) approach to the collected data. Estimations from individual low and medium‐high conversion data are compared to those found with the combined data (full conversion range data). The highest certainty in point estimates are obtained with analysis of the full conversion range data. Furthermore, the reactivity ratios determined from the combined data fall between those found with analysis of individual low and medium‐high conversion data, another corroboration of reliable data collection. Reactivity ratios determined from analysis of the combined data (r_SiA/MMA = 0.4185, r_MMA/SiA = 1.3754, r_SiA/BA = 0.8739, r_BA/SiA = 0.5736, r_BA/MMA = 0.3692, r_MMA/BA = 1.7919) are used in the recast AG model to predict cumulative terpolymer composition as a function of conversion. The experimental data and model prediction show satisfactory agreement.     

Representations for conditional expectations and applications to pricing and hedging of financial products in Lévy and jump-diffusion setting

In this article, we derive expressions for conditional expectations in terms of regular expectations without conditioning but involving some weights. For this purpose, we apply two approaches: the conditional density method and the Malliavin method. We use these expressions for the numerical estimation of the price of American options and their deltas in a Lévy and jump-diffusion setting. Several examples of applications to financial and energy markets are given including numerical examples.